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Abstract 

Observations of supernovae of type la require dark energy (some unknown exotic 'matter' of 
negative pressure) to explain their unexpected faintness. Whereas the simplest and most favoured 
candidate of dark energy, the Einsteinian cosmological constant, is about 120 orders of magnitude 
smaller than the theoretically predicted value. Motivated by this problem, a number of models of 
dynamically decaying dark energy have been proposed by considering different phenomenological 
laws or potentials of the scalar field, which are more or less ad-hoc. However, it is more advisable 
to consider the symmetry properties of spacetime rather than the ad-hoc assumptions. 

In this view, we consider a model of Robertson- Walker cosmology emerging from a Ricci sym- 
metry which provides consistently an evolving dark energy. We test the model for the recent 
supernovae la data, as well as, the ultracompact radio sources data compiled by Jackson and 
Dodgson. The model fits the data very well. 
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I. INTRODUCTION 



It is generally believed that the present expansion of the universe is accelerating due to 
the presence of some unknown cosmic 'matter' with negative pressure generally termed as 
dark energy. The simplest and the most favoured candidate of dark energy is Einstein's 
cosmological constant A, which is though plagued with horrible fine tuning problems. This 
happens due to the presence of two values of A differing from each other by some 120 orders 
of magnitudes: the magnitude of A at the beginning of inflation and the value given by 
the present-day observations. This has led a number of cosmologists to consider models of 
evolving dark energy, by either proposing different phenomenological laws or by considering 
different potentials of the scalar fields which are more or less ad-hoc. The precise mechanism 
of the evolution of dark energy, which could be required by some symmetry principle, is not 
yet known. Indeed, since the basic motivation in these models is to understand the present- 
day smallness of A, they do not provide any natural relation between the magnitude of A 
at the beginning of inflation and the present-day observational value. 

It would be worth while to investigate some symmetry principles behind the problem 
crying out for the evolution of dark energy and thus develop a more realistic and fundamental 
model for dark energy. Moreover it is always reasonable to consider symmetry properties 
of spacetime rather than considering ad-hoc assumptions. In this view, we have discovered 
a model resulting from a contracted Ricci-collineation which, apart from having 
interesting conservation properties, does provide a dynamical law for decaying A. The 
physical properties of the model have been discussed in 

HQS 

and it was found that the 

model had a credible magnitude-redshift (rn-z) relation for the observations of supernovae 
(SNe) of type la from Perlmutter et al. 4|. Since then many observations of SNe la have 
been made, quite many at higher redshifts, by the Hubble Space Telescope. It would be 
worthwhile to examine how well (or badly) the model fits the new data. 

Like the luminosity of a standard candle, the angular size of a standard measuring rod 
changes with its redshift 2 in a manner that depends upon the parameters of the model. 
Tence the Q-z relation is also proposed as a potential test for cosmological models by Hoyle 
5)]. Therefore, it would also be worthwhile to examine the O-z relation in this model for 
the dataset of Jackson and Dodgson [6|, which is a trustworthy compilation of ultracompact 
radio sources and has been already used to test different cosmological models jf], Q|. In 
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the following section we describe the model in brief 
observational relations. More details can be found in 



"or ready reference and to derive the 

a. 



II. THE MODEL 

As we are going to consider a dynamical A in Einstein's theory, it would be worthwhile to 
mention a general result which holds irrespective of the dynamics of A: the empty spacetime 
of de Sitter cannot be a solution of general relativity with a dynamical A(t) [8]. This follows 
from the divergence of the field equation [R^ — ^Rg^.j = = — g t: > A(t) /8nG . 

Obviously a solution with a dynamical A(t) is possible only if T^K ^ (and T)K. ^ 0). 
We assume that the universe is homogeneous and isotropic, represented by the Robertson- 
Walker (RW) metric, and its dynamics is given by the Einstein field equations 

S AttG 

where S and k are respectively the scale factor and curvature index appearing in the RW 
metric; and p tot = p m + Pa = p m + A(t)/8irG, p to t = Pm +Pa = Pm - Pa, with A representing 
the cosmological term. 



It is well known that collineations of Ricci tensor (Rij) have interesting s; 



mmetry prop- 
yl. The Ricci 



erties and lead to useful conservation laws in general relativity 

HQ 

collineation along a vector rji is defined by vanishing Lie-derivative of R^ along rjf C v Rij = 0. 
It has been shown that in general relativity the contracted Ricci-collineation along the 
fluid flow u % (normalized velocity 4- vector), i.e., g^C u Rij = 0, leads to the conservation of 
generalized momentum density: 

{(Ptot + 3p to tK};j = 0. (3) 

For the RW metric, the conservation law ([3]) leads to the conservation of the total active 
gravitational mass of a comoving sphere of radius S: 

(Ptot + 3p to t)5' 3 = constant = A (say) (4) 

(A detailed discussion elaborating on the meaning of (jl]) has been done in Consequences 
of the resulting models for the case A = have been discussed in [l|.) As there are 4 



unknowns p m , p m , p\ (= A(t)/8irG) and the scale factor S, the equations ([I]), (jSJ) and (jlj), 
together with the usual barotropic equation of state p m / p m = constant = w (0 < w < 1) for 
the matter source, provide a unique solution of the model. By using (jlj) in the Raychaudhuri 
equation ([T]) and integrating the resulting equation, we get 

£2 
S 2 



8nGA B 
3S 3 + S 2 ' 



(5) 



supplying the dynamics of the scale factor, where B is a constant of integration. Equations 
(J2J), (jl]) and (J3]) supply the unique dynamics of dark energy density p\ and p m as 
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(l + w) 

It must be noted that we have not assumed the conservation of the matter source which 
is usually done through the additional assumption of no interaction (minimal coupling) 
between different source fields (except for the case with a constant A which is consistent 
with the idea of minimal coupling), which though seems ad- hoc and nothing more than a 
simplifying assumption. On the contrary, interaction is more natural and is a fundamental 
principle. Let us recall that the only constraint on the source terms, which is imposed by 
Einstein's equation (through the Bianchi identities), is the conservation of the sum of all 



the energy-momentum tensors, individually they are not conserved: [R^ — \Rg 



'A 







rpij , rpij , rpij . 

i (m) + V) + + - 



implying creation or annihilation for the case A(t). 



It may be noted that the evolution of p\, as given by equation (jSJ), is a function of the 



equation of state of matter. This may 
A oc S~ 2 proposed by many authors 



3e regarded as a kind of generalization of the ansatz 
lo| . which is obtained in the present model in the 
present phase of evolution of the universe (w — 0). 

In order to study the Q-z and the m-z relations in the model, we can rewrite equations 
([MZ!) by specifying the constants A and B in terms of the cosmological parameters in the 
present phase of evolution: A = (Q - 2Q A0 )3S$H$/8>kG, B = (2fi A0 - fi + 1)H$S$, 
giving 



Qa(z) = Q A o(l + z 



2 H 

H 2 (z) 



(8) 



Q m (z) 



[(Q m0 - 20vo)(l + zf + 2fi A0 (l + z) 



H 2 



'H 2 (z) 



(9) 



where 

h(z) = H [(n m0 - 2^(1 + zf - (n m0 - 2n A0 - i)(i + z f} 1/2 . (10) 

Here fli are, as usual, the energy density of different source components of the cosmological 
fluid in units of the critical density 3H 2 /8ttG (i denoting matter (m), cosmological term (A), 
etc.). The subscript '0' denotes the value of the quantity at the present epoch. 
The angular size-redshift (Q-z) relation in the model is given by Hoyle's formula 

d(l + z) , , 

9 = v - ; radian, 11 
rS 

which relates the apparent angular diameter of the source (of redshift z located at a 
coordinate distance r) with its absolute angular size d (presumably same for all sources). 
The coordinate distance r can be calculated from the RW metric according to its curvature 
parameter k: 

siu (A So im) ■ when k = 1 
, sinh (i So im)> when k = -i. 

The present value of the scale factor So, appearing in equations ffTTl 112]) . which measures 
the present curvature of spacetime, can be calculated from 



VCOnO + fiAO-l)' (13) 

As the measured angular sizes of the radio sources in the dataset from Jackson and Dodgson 
are given in units of milli arc second (mas), we rewrite equation (TTTi) as 



Q(z; SU, dh ) = 0.0688 x dh Q sj ( ° m ° + ^ A ° ^ mas, (14) 

where d is measured in pc (par sec) and ho is the present value of the Hubble constant 
H measured in units of 100 Km s _1 Mpc _1 . We are now able to calculate the angular size 
Q(z, fli, dh ) of a radio source at a given redshift z predicted by the model (for a given set 
of the parameters fli, dho) by using equations (fTUl [T2lfT4]) . 

We also recall that the usual magnitude-redshift (m-z) relation in a homogeneous and 
isotropic model (based on the RW metric) is given by 

m{z; Qi, M) = M + 5 log{H d L (z; a)}, (15) 



where m is the apparent magnitude of a SN of redshift z located at the coordinate distance 
r, M. = M — 5 logH + 25, M is the absolute magnitude (presumably same for all SNe la), 
and g?l is the luminosity distance of the SN given by 

d L (z; Qi) = (1 + z)S r Mpc. (16) 

Now we can calculate the magnitude m(z, fij, A4) of a SN at a given redshift z predicted by 
the model (for a given set of the parameters f^, M.) by using equations ffTUl [121 [T3"l [TBI [TBI) . 
In order to fit the model to the observations, we calculate x 2 according to 

JV 

x 2 = E 

i=l 

where X° bs is the observed value of the observable, X pred (zi) is its predicted value at the 
redshift the uncertainty in the observed value X° bs and TV is the number of data 

points (or bins). (X stands for 6 and m respectively for the data of radio sources and SNe 
la.) 

III. FITTING THE MODEL TO THE RADIO SOURCES DATA 

We use the sample of 256 ultracompact radio sources compiled by Jackson and Dodgson . 
This sample of 256 radio sources with z in the range 0.5 to 3.8 was selected by them from a 
bigger sample of 337 ultracompact radio sources originally compiled by Gurvits These 
sources, of angular sizes of the order of a few milliarcseconds (ultracompact), were measured 
by the very long-baseline interferometry. The points of the sample of Jackson and Dodgson 
are short-lived quasars deeply embedded inside the galactic nuclei, which are expected to be 
free from evolution on a cosmological time scale and thus comprise a set of standard rods, 
at least in a statistical sense. Jackson and Dodgson binned their sample into 16 redshift 
bins, each bin containing 16 sources. We fit the present model to this sample by calculating 
X 2 according to ffTTj) and minimize it with respect to the free parameters f2 m0 , Q\o and dh Q . 
The global minimum is obtained for the values (with the constraint f2 m0 > 0) 

tt m0 = 0.68, tt A0 = 0.45, dh = 6.50 with X 2 = 13.00 

at 13 degrees of freedom (DoF). This represents a very good fit, with the goodness-of-fit 
probability Q = 44.8% (see the Appendix for an explanation of Q). We further note that the 
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FIG. 1: Some best-fitting models to the ultracompact radio sources data from Jackson & Dodgson 
are shown. The solid curve corresponds to the best-fitting dark energy model pa oc S~ 2 with 
^mo = 1 ~~ ^Ao = 0.62, the dashed curve corresponds to the best-fitting standard ACDM model 
with f2 m o = 1 — ^ao = 0.21 and the dotted curve corresponds to the best-fitting Einstein-de Sitter 
model (f2 m o = 1 — ^ao = 1)- 

solution for the minimum \ 2 is very degenerate in the parameter space and the parameters 
wander near the global minimum of x 2 m almost a flat valley of some complicated topology. 
For example, the best-fitting flat model (ft\o = 1 — fi m o) is obtained as 

Q m0 = 0.62 ± 0.03, dh = 6.46 ± 0.31, x 2 /DoF= 13.09/14, Q = 51.9%, 

which represents a slightly better fit (as the number of DoF is increased). It may be 
mentioned that the value f2 m o = 0.62 estimated from these data is in good agreement 
with the results estimated (in the following section) from the recent observations of SNe 
of type la. It should be noted that the best-fitting solution Q m0 = 0.62 gives a mildly 
accelerating expansion of the universe at the present epoch: the deceleration parameter 
g = Q m o/2 — f^Ao = —0.07. In order to compare, we find that the best-fitting concordance 
model (flat standard ACDM model with a constant A) to this dataset is obtained as 

n m0 = 0.21 ± 0.08, dh = 7.25 ± 0.55, x 2 /DoF= 16.03/14, Q = 31.1%, 
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which also represents a good fit. One may note that the value f2 m o = 0.62 obtained for 
the model pa oc S~ 2 is higher than Q m o = 0.21 obtained for the standard ACDM model. 
However, one should note that the other precise observations, which one would expect to 
be consistent with any model, are the measurements of the temperature anisotropy of CMB 



made by the WMAP experiments [12j, whose only apparent prediction is f2 m0 + fiAo — 1 [13 ] . 
For this reason, and also motivated by theoretical considerations required by inflation and 
flatness problem, we assume spatial flatness henceforth. For curiosity, we test the Einstein- 
de Sitter model (fl m o = 1, A = 0) against the radio sources data: the best fitting solution is 
obtained as 

dh = 4.91 ± 0.11, x 2 /DoF= 28.82/15, Q = 1.69%, 

which can be rejected by the data only at 98.3% confidence level. In Figure 1, we have 
shown some models obtained from our fitting procedure and compared them with the data. 

IV. FITTING THE MODEL TO THE RECENT SNE IA DATA 

It was already shown that the model explained the data of SNe of type la from Perlmutter et 
al. successfully 4j. Since then many SNe of type la at higher redshifts have been discovered. 
Extending our earlier work further, we now examine how the model fits the updated gold 
sample of Riess et al. [ijj]. In addition to having previously discovered SNe la, this sample 
(of 182 SNe la in total) contains 23 SNe la at z > 1 recently discovered by the Hubble Space 
Telescope and is claimed to have a high-confidence-quality of spectroscopic and photometric 
record for the individual SNe. 

We find that the present model has an excellent fit to this data, comparable with the 
concordance model. The best-fitting concordance model is obtained as 

tt m0 = 0.34 ± 0.04, M = 43.40 ± 0.03, x 2 /DoF= 158.75/180, Q = 87.1%, 
an excellent fit indeed! The model A oc S~ 2 provides a similar fit: 

tt m0 = 0.58 ± 0.02, M = 43.46 ± 0.03, x 2 /DoF= 167.90/180, Q = 73.2%. 
The Einstein-de Sitter model does not fit the data well: 
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FIG. 2: Modified Hubble diagram of new 'gold sample' of SNe la [1J] minus a fiducial model 
(f2 m o = Oao = 0): The relative magnitude(Am = m — wifiducial) is plotted for some best-fitting 
models. The solid curve corresponds to the dark energy model p\ oc 5~ 2 with O m o = 1— ^ao = 0.58, 
the dashed curve corresponds to the standard concordance model with f2 m o = 1 — = 0.34 and 
the dotted curve corresponds to the Einstein-de Sitter model. 

M = 43.72 ± 0.02, X 2 /DoF= 283.40/181, Q = 1.7 x 10" 4 %. 

These models have been shown in Figure 2. In order to have a visual comparison of the fits 
of different models to the actual data points, we magnify their differences by plotting the 
relative magnitude with respect to a fiducial model f2 m0 = ^ao = (which also has a good 
fit: x 2 /DoF = 174.29/181, Q = 62.6%). 

V. CONCLUSIONS 

In order to test the consistency of the cosmological models with observations as well as to 
estimate the different cosmological parameters, data on SNe of type la and radio sources 
have been used by several authors. We use the recent gold sample of SNe la from Riess et 
al. and the sample of 256 ultracompact radio sources (of angular sizes of the order of a few 
milliarcseconds) compiled by Jackson and Dodgson to test a model of dark energy which 
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results consistently from the contracted Ricci-collineation along the fluid flow vector. We 
find that the model has excellent fits to both the data sets and the estimated parameters 
are also in good agreement. 
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APPENDIX 

Though there is not a clearly defined value of y 2 /DoF for an acceptable fit, however it 
is obvious from equation (jTTj) that if the model represents the data correctly, the difference 
between the predicted angular size /magnitude and the observed one at each data point 
should be roughly the same size as the measurement uncertainties and each data point will 
contribute to x 2 roughly one, giving the sum roughly equal to the number of data points 
TV (more correctly N— number of fitted parameters = number of degrees of freedom 'DoF'). 
This is regarded as a 'rule of thumb' for a moderately good fit. If x 2 is large, the fit is 
bad. However we must quantify our judgment and decision about the goodness- of- fit, in the 
absence of which, the estimated parameters of the model (and their estimated uncertainties) 
have no meaning at all. An independent assessment of the goodness-of-fit of the data to the 
model is given in terms of the x 2 -P r °bability: if the fitted model provides a typical value of 
X 2 as x at n DoF, this probability is given by 



Roughly speaking, it measures the probability that the model does describe the data genuinely 
and any discrepancies are mere fluctuations which could have arisen by chance. To be more 
precise, Q(x, n) gives the probability that a model which does fit the data at n DoF, would 
give a value of x 2 as large or larger than x. If Q is very small, the apparent discrepancies 
are unlikely to be chance fluctuations and the model is ruled out. For example, if we get 
a x 2 — 20 at 5 DoF for some model, then the hypothesis that the model describes the data 
genuinely is unlikely, as the probability Q(20, 5) = 0.0012 is very small. It may however 
be noted that the x 2_ P r °bability strictly holds only when the models are linear in their 
parameters and the measurement errors are normally distributed. It is though common, 
and usually not too wrong, to assume that the x 2 - distribution holds even for models which 
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are not strictly linear in their parameters, and for this reason, the models with a probability 
as low as Q > 0.001 are usually deemed acceptable fis| . 
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